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The integrals in Gradshteyn and Ryzhik. 
Part 16: Complete elliptic integrals 

f^ , Stefan Boettner and Victor H. Moll 

(N 

^^' Abstract. The tabic of Gradshteyn and Ryzhik contains many entries that are 

^ related to elliptic integrals. We present a systematic derivation of some of them. 

Si 

1. Introduction 

Elliptic integrals were at the center of analysis at the end of 19*''-century. The 
complete elliptic integral of the first kind defined by 

dx 



< 



a 



V(f-a;2)(l-fc2a;2) 

is a function of the so-called modulus k^. The corresponding complete elliptic 
integral of the second kind defined by 



(1.2) EW:.^'y™dx. 

^~^ . The total collection of complete elliptic integrals contains one more, the so-called 

^2^ I complete elliptic integral of the third kind defined by 

CN ' f^ dx 

(f.3) n(n,fc):=/ =. 

^— N ■ The complementary integrals are defined by 

(f .4) K'(fc) - K(fc') 



where k' — \/l--W is the so-called complementary modulus. 
^^ . The change of variables x = sin t yields the trigonometric versions 

r^/2 u riT/2 

(f.5) K(fc) = / - and F,(^) = / Vl-k^sin^tdt, 

Jo \/l - k"^ sin^ t Jo 
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2 S. BOETTNER AND V. MOLL 

with a similar expression for II(n, k). 

In general, an elliptic integral is one of the form 

(1.6) ^^ /^P(^^ 

Ja V 

where y^ is a cubic or quartic polynomial in x. The integral is called complete if a 
and h are roots of y = 0. It is clear that K(fc) is elliptic. The same is true for E(fc), 
written in the form 

, , :^ (l~k^x^)dx 

(1.7) E(A)- ' ^ ' 



'O y/il-X^){l-k^X^) 

2. Some examples 

In this section we offer some evaluations from 4^ that follow directly from the 
definitions. Some special values are offered first. The evaluations of these integrals is 
facilitated by Legendre's relation 

(2.1) K(fc)E'(fc) + K'(fc)E(fc) - K(fc)K'(fc) = |. 
The reader will find this identity as Exercise 4 in section 2.4 of [7]. 
Example 2.1. 

(2.2) K(V^) = ^=r^(- 

The proof is direct. The integral is 

^ ' ^ ' Jo VT^^ 4io ^ ^ 4r(3/4) 

The result now follows from the symmetry rule 

(2.4) r(a)r(l - a) = ^^ 

sm 7ra 

for the gamma function and the special value r(l/2) = y/ir. This example appear as 
entry 3.166.16 in W. Entry 3.166.18 states that 

r2 5) f' ^J^ = J_Y^(- 

^ ' ^ Jo VT^^ V2^ V4 

The proof consists of a reduction to a special value of the beta function. The change 
of variables t — x"^ gives 

This integral is ■jB (|, i). The simplified result is obtained as above. 
Formula p.7p with k = \/— 1 shows that 

(2.7) E{V^)^J^ -7=^dx. 



VT 
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The values given above show that 
(2.8) E(^/^) = ^ 



rM-)+4rM- 

4\/27r 
Example 2.2. 

This appears as entry 8.129.1 in 4 . This value conies from the previous example 
and the identity 

(2.10) K(V^fc) = -jl=K ^ ^ 



VT+fc2 V\/l + A:2 



with fc = 1. The identity (|2.10p follows by the change of variables x i-> x/ ^1 + fc2(l — x^) 
in the left-hand side integral. 

The values of the modulus k for which K'/K is the square root of an integer are 
of considerable interest. These are called the singular values. The previous example 
shows that 1/V2 is the simplest of them: in this case 

fa)-' 

A list of the first few values fc^ for which 

(2.12) ^(fc,) = Vf 

is given in [3] and it starts with 

fe = V2 - 1, fcs = ^^"^\ fc4 = 3 - 2V2, fc5 = ^ NVI-I - \/3-%/5 

3. An elementary transformation 

Elementary manipulations can be emploted to evaluate certain entries in [4] . For 
instance, direct integration by parts on the integrals defining the functions K and E 
produces 

/■ X arcsin x , 1 / vr , , , , , 

and 

/„ ^N f a; arcsin x , 1 /^,,n tt , ,\ 

This last evaluation appears as entry 4.522.4 in |4j. 

On the other hand, several entries in [1] may be evaluated also by integration by 
parts choosing the inverse trigonometric term to be differentiated. Such procedure 
gives 
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_ „, f X arccos x da; 1 /tt „,,,\ 

that appears as entry 4.522.5, 

,„ ,, f xarcsinxdx 1 /tt „,,, 

that appears as entry 4.522.6, and finally 4.522.7: 

,„ ^, /"^ a;arccos xdx 1 / tt , , .^,,, 

In this section we derive a different type of elementary transformation for integrals 
and use it to obtain the value of some elliptic integrals appearing in [4j . 

Lemma 3.1. Let f be an odd periodic function of period a. Then 

(3.0) ri^^.iri^,,^ 



jQ -v u jQ tan — 

Proof. The result follows by splitting the integral as 

a/2 






X + ka (fc + l)a — X 



dx 



k=0 

and using the partial fraction decomposition 

.o 7N . ^b 46 ^ 1 

(^■^) ^^'^y-V^^ (2^-1)2-52' 

given as entry 1.421.1 in j4]. D 

Corollary 3.1. Let f be an even function with period a. Then 

(3.8) r li^sin — dx^- r f{x)dx. 
Jo X a a Jq 

In particular, for a = n, 

/■oo y^^-j ^17/2 

(3.9) / sin a; do; = / f{x)dx. 

Jq X Jq 

Proof. Apply the lemma to the function f{x) sin ^^ which is odd and it has 

period 2a. The result follows from the half-angle formula 

X sin a; 

3.10 tan- = 

2 1 + cos a; 

and the value 

/■" nx 

(3.11) / f(x)cos — dx^O. 

Jo a 
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n 

A similar results holds for odd functions. These appear as entry 3.033 in [4|. 
Corollary 3.2. Let f be an odd function with period a. Then 

(3.12) / sin — dx — — I f{x)cos — dx. 

Jo X a a Jf^ a 

In particular, for a = it, 

f°° fix) /■'^^^ 

(3.13) / si'nxdx= / f{x)cosxdx. 

Jo X Jq 

Example 3.1. The function f{x) = 1 and a = vr in Corollary 13. II gives the classical 
integral 



(3.14) r= 

Jo 



oo 

sma; , n 
dx — —. 



iQ X 2 

This is entry 3.721.1 in |^. The reader will find in [5l [6] a couple of articles by G. 
H. Hardy with an evaluation of the many proofs of this identity. These papers are 
available in volume 5 of his Complete Works. 

Example 3.2. Entry 3.842.3 of [4] consists of four evaluations, the first of which 

(3.15) r -'"^^ f^K(.). 

"'0 V 1 - fc^sm X X 



This follows from Corollarv 13.11 bv choosing a — n and f{x) = 1/vl — k^sin^ x. A 
different proof of this evaluation is offered in Section [5] below. 

Example 3.3. A second integral appearing in 3.842.3 is 

/„ „x f°° sina; dx .,^,,, 

(3.16) / ,^ ,, , — ^K(fc). 

Jo V 1 - k^ cos2 X X 

also follows from Corollarv 13. II This is also true for entry 3.841.1 

/>00 7 

(3.17) / sin a; Vl - fc2 sin^ a; — = E(fc) 

Jo X 

and its companion entry 3.841.2 

(3.18) / sinxVl-fc^cos^x— = E(fc). 
Jo X 

Example 3.4. The elementary method introduced here may be used to evaluate all 

integrals of the type 

f°° sin" X cos" X dx 
(3-19) /m,«(fc) := / , = — 

"'0 V 1 — fc sin X ^ 

and the companion family 

,„ „„, ^ ,,, f°° sin" X cos™ a; dx 

(3.20) J„,,„(fc) := / , ,, , — . 

Jo V 1 - k^ cos2 X X 

All entries in the sections 3.844 and 3.846 match one of these forms. 
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Example 3.5. Many other evaluations can be produced by this method. For instance, 

,„ „ , f°° sinx log(l - P sin^ x) dx r^'^ logfl - P sin^ x) , 

(3.21) / = = / — ax. 

Jo V 1 - fc^ sin^ X ^ Jo \J\-'k? sin^ x 

The integral on the left appears as entry AAZ1.\ and the one on the right is entry 
4.414.1 in [4]. A proof of the identity 

(3.22) r^Ml=£l£^ci. = K(fc)ln.', 
Jo Vl — k"^ sin^ X 

is given in Example 17.21 

4. Some principal value integrals 

The method described above can be employed to evaluate some entries of |4] 
provided the integrals are interpreted as Cauchy principal values. 

Example 4.1. The first example is 

(4,1) r /, T. f-'^w- 

that appears as one of the four entries in 3.842.3 of [4]. 

Let I\ (fc) denote the integral and introduce the notation 

(4-2) m = -^^- 

V 1 — fc sin X 
Then j(x) is odd and it has period tt. The principal value of the integral is given by 

(4.3) h{k) = lim y I /"^' ' JMdx+ r -l^ dx] . 

The substitution y = vr — a; in the second integral above produces 



X x — (j + l)7r 





/i(fc) 



lim , I ' / o .o 9 

The series corresponds to the partial fraction expansion of the cotangent function. 
This completes the evaluation of (j4.ip . The reader will note that this proof is very 
similar to that of Lemma 13.11 
The value 

(4-4) / — = K(fc), 

Jo V 1 - fc COs2 X X 

that also appears in 3.842.3 is established using the same type of argument. This 
completes the evaluation of the integrals in that entry of [4] . 
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Example 4.2. Entry 3.841.3 of [4] 

(4.5) / tan a; V 1 - fc^ ^^^ x—^ E(fc) 

h X 

and its companion 3.841.4 

roc ^^ 

(4.6) / tanxVl - fc^ cos^ cc — = E(A:) 
Jo X 

can be established by the method described in the previous example. 



5. The hypergeometric connection 

The identites among elliptic integrals often make use of the series representations 

1 



(5.1) 


K{k)^-2F, 


and 

(5.2) 


E(fc) = |2Fi 



1 

2 2-1,2 



1 1 

'2 2. /,2 



9 ^ 



i=o 



IE' 



J'- 



i\ ' 






i=o 



where 2-F1 is the classical hypergeometric function 

a b 



(5.3) 



iFi 



X 



V (°)j (^)j ^.7 



and 

(5.4) (a), = a(a + l)(a + 2) • • • (a + J - 1), 

is the Pochhammer symbol. The value (a)o = 1 is adopted. 

6. Evaluation by series expansions 

In this section we describe a method to evaluate many of the elliptic integrals 
appearing in [4]. 



Example 6.1. The first example is entry 3.842.3 

sin a; 



(6.1) 



\Jl - k? sin^ X 



— = K(fc) 

X 



that has been evaluated in Section [S] 
Define 



(6.2) 



/i(fc') := 



sma; 



dx 







Vl-fc^ 



sm X 



To evaluate the integral, let m — k^ and expand the integrand in power series using 

sin a; / 1\ sin •' x 



dm 



vv 



X 



-(1 — msin x) 



-1/2-j 
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Therefore, 

The remaining integral is entry 3.821.7 in [J: 

The value of the integral (j6.ll) now follows from the series representation of K(A:) given 
in (ED. 



Proof of dm. Start with 

(6.6) sin^J+i X ^ 2-2J ^{-ly-'^ (^^ + ^"j sin(2j -2y+ l)x 



and the integral in Example 13.11 in the form 

f°° sin ax tt 

for a > 0, to obtain 

(6-8) X ^~'' = ^^S"'^ 

It follows that 



°° a\ {-lym^ 4^, .^.^2j + l 



The result now follows from the next lemma. 
Lemma 6.1. Let j, fc e N. Then 

(6.10) j:{^iy{^'^') = {-if(^'^ 

u=0 V ^ / 

Proof. The proof is by induction on k. The case fc = is clear. The induction 
hypothesis is used to produce 

(0.11) i:(-irfr)^(-i.-G^o.(-i)'e^r' 

and an elementary calculation reduces this to (— l)'^(^j^). This completes the proof of 
(ESD- □ 

Second proof of (|6.5p : apply the identity p.l3p to the function f{x) ~ sin ^ x to 
obtain 

(6.12) / — -dx= sin^^xdx. 

Jo ^ Jq 
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This last integral is the classical Wallis' formula given by 
(6.13) / sm''^xdx ^^'' 



JO 2 j\ 

The reader will find in ^1^ information about this formula. 

Example 6.2. Entry 3.841.1 in [i] 

/•OO 7 

(6.14) / sin a; Vl - fc2 sin^ 2: — = E(fc) 

Jo 2; 

is established by the same method employed above. The proof starts with the expan- 
sion of the integrand using 

^^irN f d, V sina; r r^— . .-. sin^-'+"^ a; , . . 2 ni/2-i 

(6.15) -;— VI -TOsm'^x = (-i) . (1 - msm^x)^/"^ ^ 



and then identify the result with (157 

Example 6.3. Entry 3.842.4 in [4] states that 

_ „, ^ ,,, (^''^ a; sina; cos a; , Trfc' E{k\ 

(6-16) /2(A;):=/ -====da: = -— + ^. 

"'0 V 1 - fc^^sm a; ^k, /t 



The parameter k' is the complementary modulus fc' = \/l — fc^. 
Write m = k^ and expand the integrand in series using 

, d \'' a: sin a; cos a; /l\ xsin-*^ a: cos a; 

(6.17) 



^dm, J yl — 771 sin x V / j vl^^T-sin x 
Therefore 

(6.18) hi'm) = 7 I - I — ;- / ajsin^^"*""^ a; cosxdx. 
Integration by parts gives 

(6.19) [''-'^''^' - -«-^- - i(-^ - i(-^^ (^- + i' ^) ' 
where 

(6.20) B(w, u) = / r-i(l - t)'"-'^ dt = 2 sin^"-! (^ cos^^-V^'^, 

Jo Jo 

is the classical beta function. It follows that 

,6.21, .,,„, = -g(-)^^_^--g(-)^,,(, + -.-)^_^. 

The two series are now treated separately. 
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The first sum is computed by tlie binomial tlieorem 



.22) 



OO , N 



3=0 



as 



(6.23) 



4 2^ \ 9 



4^V2;,(J + 1)! 2(l + yr^7^) 2to 



(1 - Vi^^). 



The second sum is 



4-^ V 2 / . 



2' 2; (j + 1)! 



■E 



ra 



4 ,t^V2y,(j + l)!r(j + 2) 






8^V2yA2A (j + 1)! 



8 ^ 

-y - 



2m ^^ \ 2 , ., , 



.j+i 



2A+i(j + l)! 



2m 



^Fi ( 2 ^ 2.^1-1 



The hypergeometric representation (j5.2p and (j6.2ip give 

/2(m) 



(6.24) 

as claimed 



7rVl-TO E(fc) 



2m 



7. A small correction to a formula in Gradshteyn and Ryzhik 

In this section we present the evaluation of some elliptic integrals in [4] . In par- 
ticular, a small error in formula 4.395.1 is corrected. 



Proposition 7.1. Let k' = y/l — fc^ ^g the complementary m,odulus. Then 

In a; da; 1 , 



(7.1 



^{l + x^){k'^+x^) 2 
Proof. Let m = k'^ and use 



K(fc) Ink'. 



(7.2) 1:^'' 



Inx 



(-IF 7T 



In a; 



2Jj ^(l + a;2)(m + x2)J+i/2 



,c?TO/ ^(1 +a:2)(m + a;2) 
to expand the integrand around m, = — 1. It follows that 

In a; da; v-^ (— l)-* /in f°° In a; da; 



(7.3) 



, == V (~^F (1) / ^"""" fm-1)^ 

V(l + ^')(^''+^') ,t^o ^^ '^0 (l + x2)^+i^ '- 
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is last integral is given by 




/■°° \nxdx 


1 r°° Inxdx 


Jo (l + x2)^+i 


4 Jo yi(l + a:)J+i 




= 4da^("'^ " + ^) 




-L„/i . i\r,/i\ 



= 1/2 



B{hj + '.)[^{'2)-^{j + '.)] 



_ n n 

■' ^ ■' J i=0 
Therefore, the left-hand side of (17.11) satisfies 



'^ 2i + l 



(7.4) 



In a; dx 



oo (lY J-i 



V(l + '^')(^'"+2^') 



^ 2^ j!2 Z^2^ + l^ ^ 



The series expansion for the complete elliptic integral now shows that the right-hand 
side of (|7.ip is given by 



1 



In ?7iK(\/l — rn) 



E 

i=i 



(1 - my 



oo n)2 



■E 

J=0 



A — n -^ 



.=0 J- 
(l-m^. 



The result follows from the identity established in the next lemma. 
Lemma 7.1. Let j G N. Define 



n 



(7.5) 
Then 



i-i 



r\2 



i-1 



1=0 -^ i=0 

Proof. The relations 
(7.7) (-X). = (-l)n.x - k + 1), and (i)„_, (i - n), = (-1)'= (i)„ 

can be used to rewrite the left-hand side as 



i-i n\2 



E^ 



(« 1 



j-i Ci^ 



i!^ J — i 



E 

fc=0 



1 



.2yj-fc-i 

{j -k~lP k + l 
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Thus the assertion of the lemma is equivalent to 



(7.8) 



1 



E \ J)k+1 

k=0 \2 ■'Jk+l 



J-1 



T=i: 



2k + 1 



'k+l fc=0 

Next apply the fact that {x)k+i = a;(x + l)fc to obtain 



1 



E V J)k+1 

(l-n)^ k + 

fc=0 \2 Jlk+l 



J 



3-1 

E 



(i-j?fe 1 






•1 a\^^' 



The right-hand side is a balanced /^F^ series and it can be transformed using 



4^3 



y 



u V 



■ 4F3 



W 
U — X 



u — y z 

1 — V + z — m 1 — w + z — m u 



See [2], page 56. Now let y = z ~ 1, x ^ 1 — j, m = j — 1, u = v ^ ^ — j and w ~ 2. 
It follows that 

"1 



4-F3 



3 _ • 3 ■ - ^ J^ 

2 J 



2 J 
(l-jVl(l)i-l 

(|-jVi(2),-i 



4-F3 



1 



-^•+1 



-J + 



1-J 



1-3 



The last hypergeometric terms is now simplified 
_ (2j-l) 

k 



2^--ig(^). H + i). 



J 



k=o V2;; 



-3 + 



2Jk 



3 




(2j- 


1)2 


2j' 


^ 


(2.7 - 


1)2 



2 J-l 

E 

fc=0 

E 

fc=0 

i-1 



1 



{2k + l){2j - 1 - 2k) 



1 



1 



2fc + 1 2j -l-2k 



1 
^2fc + l' 



fc=0 



as claimed. 



D 



An automatic proof. The result of Lemma 17.11 also admits an automatic proof as 
described in ^8". Define the functions F{i,j) and G{i,j), respectively, as 



(7.9) 



Fi^,J) 






1 



andG(*,j) = -^il|^ 



2Jj + l *• 
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The stated result is equivalent to the identity a{j) — b{j), where 

(7.10) aij) = ^^(*,j) and 6(j) = ^ y--. 

Zeilberger algorithm finds the non-homogeneous recurrence 

(7.11) F{t + l,j) - F{i,j) - G(^ + l,j) - G{i,j). 

Summing this for i from to j — 1 and using the telescoping of the right-hand side, 
produces 

j—i j—i j—i j—i 

J2f{z,j + 1)-J2f{z,j) = ^G(z + 1,j)-EG(*'^') 

i=0 j=0 j=0 j=0 

= GO-,j)-G(0,j) 

"(2j + 1)2- 
Now observe that 

4(J + 1)^ ^-^ 



a(j + l)-a(j) = ^:^^+^F(*,j + 1)-EF(*,j) 
4(j + 1)2 Af 



(2j + 1)2 (2j + 1)2 
4 



2j + l 

The sequence b{j) satisfies the same recurrence. Therefore a{j) — b{j) is a constant. 
Since a(l) = 6(1) = 4 this constant vanishes. This establishes the result. 

The next result corrects entry 4.395.1 in [4]. 

Corollary 7.1. The value 

/■°° Intan6ld6l 1, ,,^,,,, 

(7.12) / — = — lnk'K{k) 

Jo V 1 - fc2 sin2 6* 2 

Proof. Let x i-^ tang in (fTl} . D 

Example 7.1. Entry 4.242.1 states 



(7.13) r^J^^i^^lKh^^IH^ ln«6. 

70 V(a2+a;2)(52 + a;2) 2a \^ a ^ 

Formula (j7.ip corresponds to the special case a — 1. The change of variables x ^ at 
produces 

Inxdx 1 f°° Intdt Ina f°° dt 

-- + - 



V(a2+a;2)(62 + 2;2) Mo ^(1 + ^')(c' + *') ^ ^ ^(1 + t2)(l + c2i2) 
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with c = h/a. The first integral is evaluated using (|7.ip and let t = tan (^ to see that 
the second integral is K(Vl — c^)- This establishes the result. 

Example 7.2. The techniques illustrated here are now employed to prove entry 
4.414.1 in [J: 

r^^ ln(l - A:^ sin^ x) , ^,,,,, ,, 
J a v 1 — fc2 sin^ a; 

Let m = fc^ and observe that 



d aj + (ijh\{l — msiv? x) . 2, a^+i + /3j+i ln(l — msin^ a;) . 2J+2 

dm (1 - msin^a;)J+i/2 (1 - msin^ x)-''+3/2 



where the parameters aj, Pj satisfy 

(7.16) a,+i = {j + i)a - j - /3, and /3j+i = (j + i)/?,. 

Now choose ag = and /3o = 1 to obtain 

d ^■^ ln(l — rnsin x) a^ + /3j ln(l — ?7isin x) . 2,- 



Vrf"^/ Vl-msin^a: (1 - msin^ x)J"+i/2 

Expand the integrand of (j7.14p around m — and use 

r°° (-) 

(7.18) / sm^'xdx = ^-'' 



2 jP 
and the expressions 

'"' 2 



to see that 



r/' ln(l - P sin2 x) , ^ (^) , /^ 1 



j=0 •' \i=0 



The result now follows from the evaluation given in the proof of Proposition [73] 
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